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Determination of the accurate astrophysical S factor of ^^C(q:, 7)^®0 reaction has been regarded 
as a holy grail of nuclear astrophysics for decades. In current stellar models, a knowledge of that 
value to better than 10% is desirable. Due to the practical issues, tremendous experimental and 
theoretical efforts over nearly 50 years are not able to reach this goal, and the published values 
contradicted with each other strongly and their uncertainties are 2 times larger than the required 
precision. To this end we have developed a Reduced R-matrix Theory, based on the classical R- 
matrix theory of Lane and Thomas, which treats primary transitions to ground state and four 
bound states as the independent reaction channels in the channel spin representation. With the 
coordination of covariance statistics and error propagation theory, a global fitting for almost all 
available experimental data of system has been multi-lteratively analyzed by our powerful code. 

A reliable, accurate and self-consistent astrophysical S factor of ^^C(a,7)^®0 was obtained with a 
recommended value Stot (300) = 162.7 ± 7.3 keV b (4.5%) which could meet the required precision. 

PACS numbers: 25.55.-e,25.40.Lw,26.20.Fj,21.10.-k 


I. INTRODUCTION 


During core-He burning 3a and ^^C(a,7)^®0 reactions 
compete to determine the helium burning time scale, to¬ 
gether with the convection mechanism, and the relative 
abundances of oxygen and carbon prior to core-C burn¬ 
ing. W. A. Fowler, a Nobel laureate in physics in 1983, 
definitely held a view that the abundance ratio of to 
and the solar neutrino problem are serious difficul¬ 
ties in the most basic concepts of nuclear astrophysics [ij . 
The work of many investigators has resulted in a knowl¬ 
edge of the predicted reaction rate of 3a process within 
about 10% [2| accuracy at the usual helium-burning tem¬ 
peratures. Unfortunately, the same of the ^^C(a,7)^®0 
reaction is much less well determined, and it is highly de¬ 
sirable to know this rate with comparable accuracy (3| to 
that of the 3a process in order to provide adequate con¬ 
straints on stellar evolution and the synthesis of elements, 
e.g. the yield of the neutrino-process isotopes ^Lij_/^B, 
i 9F, and ^®9Ta in core-collapse supernovae [1, 0, 

and the production of the important radioactive nuclei 

26Al, 44rpj^ 60pg |]_ 

The most direct and trustworthy way to obtain the 
astrophysical S factor of the ^9c(a,7)^®0 reaction is 
to measure the cross section for that reaction till as 
low energy as possible, and to extrapolate to energies 
of astrophysical interest. The astrophysical S factor for 
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^9C(a, 7)^90 reaction reaction is given by 

S{E) = a{E) ■ E ■ exp(27r77), (1) 

where rj = ZaZce^/hv is the Sommerfeld parameter 
for the interaction particles. To investigate the spe¬ 
cific role of the nucleus for the S factor, a wealth 
of experimental data have been accumulated over the 
past few decades, including the precise measurements 
of total cross-section of ^9c(a, 7)^90 [M , 7 -rw an¬ 
gular distributions of ground state transition (lll - [^ . 
cascade transitions [l^, HI, Hbd-delayed a spectra 
for , transfer reaction j29l - [^ , elastic scatter¬ 

ing ^9c(a,a)^9c (3l - [^ and additional particle reaction 


TABLE I: Comparison of the astrophysical S factor (keV b) 
at 0.3 MeV obtained in various fits, including this work, for 
the ElO, E20, cascade transitions components, as well as the 
total. 


Reference 

Stot 

Seio 

Se20 

Seas 

This work 

162.7±7.3 

98.0T7.0 

56.0T4.1 

8.7±1.8 

Schiirmann 

161T19 

83.4 

73.4 

4.4 

Oulebsir 

175T63 

100T28 

50T19 

— 

Sayre 

— 

— 

62/^ 

— 

Tang [2^ 

— 

84±21 

— 

— 

Matei [^j 

— 

— 

— 

S6.92=7.0±1.6 

Matei [24] 

— 

— 

— 

5'6.o5=25/?® 

Hammer [^] 

162T39 

77T17 

80T20 

4±4 

Tischhauser [3^ 

149±29 

80T20 

c:o+13 

0O_i3 

16T16 

Kunz [lA] 

165±50 

76T20 

85T30 

4±4 

Brune [29] 

159 

101T17 

42/^^ 

16 

Ouellet [IR] 

120T40 

79T16 

36±6 

— 

Buchmann \4^ 

165±75 

79T21 

70T70 

16T16 
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pathways ^^C(a, q;i)^^C and ^^C(a,p)^^N [s^l^ 

at high energies. However, complexity of reaction mech¬ 
anisms of ^‘^C+a makes it extremely difficult to deter¬ 
mine the S factor, despite five decades of experimen¬ 
tal investigations yet, the desired accuracy and precision 
associated with the ^^C(a,7)^®0 reaction continues to 
be an obstacle (See the Table IB . Recent two collabora¬ 
tions [4^, and our team have pursued comple¬ 
mentary approaches to obtain the inverse ^®0(7, a)^^C 
reaction to energies lower than the currently achieved, 
which could offer significant advantages over traditional 
approaches. But the expected outputs of these proposals 
in low energy are far from what’s required in the stellar 
models. R-matrix analysis is the most effective method 
for the fitting and extrapolation of existing data of 
system, which are main content of this work. 

The classical R-Matrix theory of Lane and Thomas [ 1 ^ 
deduced the standard R-matrix formulae to describe two 
body nuclear reaction. However, these formulae were 
thought not justified for 7 radiative capture, because of 
the possibilities of particle production and annihilation, 
and it is hard to select suitable channel radius for the 
long-range electromagnetic transition. After that, an¬ 
other paper of Lane | 49 l | expanded the collision matrix to 
radiative capture for a sum of three parts, viz. an internal 
resonant, an external resonant and a non-resonant part 
corresponding to the channel integral from hard sphere 
scattering. Based on this conclusion, the angle-integrated 
cross section formulae are derived from perturbation the¬ 
ory in Refs. [laiiii, and the adjustable parameter of pho¬ 
ton reduced-width amplitude can be split into internal 
and asymptotic channel contributions in practical appli¬ 
cations [^. Recently, a vital progress in R-matrix 
code, AZURE was presented in Refs, and it al¬ 

lows simultaneous analysis of the integrated and differ¬ 
ential data for the electromagnetic transition. 

In the most widely used theory [s^, since the S fac¬ 
tor of the El and E 2 multipoles has different energy 
dependence, one must have an independent and precise 
information on each multipole cross section for an extrap¬ 
olation to 0.3 MeV. So the secondary data of El and E 2 
multipoles were in general used for the R- matrix analy¬ 
sis in the previous publications. The primary data most 
often consist of angular distributions measured at many 
discrete energies. Each primary distribution was then an¬ 
alyzed independently in terms of the appropriate set of 
Legendre polynomials 0, a non R-matrix analysis ne¬ 
glecting the energy- and angle-dependence, to yield the 
secondary data, aEio and (Jeiq/<Je20 at this discrete en¬ 
ergy. And these data and their error values, being derived 
quantities, are no longer proportional to the experimen¬ 
tally measured quantities, i.e., the angular distribution 
yields, and thus lead to complications and discrepancy 
for the extrapolation of the ^^C(a,7)^®0 S factor | 43 |. 

A R-matrix code for a nuclear system is in princi¬ 
ple an exact model as long as a complete set of quan¬ 
tum states of a nuclear system and all the corresponding 
experimental data can be accurately described simulta¬ 


neously (called Global fitting). Any inconformity that 
does not meet the principle will induce inestimable un¬ 
certainty, because each channel, each level and each data 
are intimately correlated and strong interference in the 
nuclear system. By now the R-matrix approach men¬ 
tioned above [S^, has got popular application with the 
procedure, but their results have very larger difference 
(See Table HI). And no one used it to do a global fitting 
for system in the astrophysical energy yet. Eor these 
reasons, based on the theory of Lane and Thomas [i^, we 
develop a Reduced R-matrix Theory to make the global 
fitting for the special problem to search for the S factor 
of i2c(a, 7)160. 

Section H summarizes the construction of the Reduced 
R-matrix Theory and general aspects of the R-matrix 
approach in the global analysis. Section HI presents the 
construction of reaction channel, evaluation and fits of 
experimental data. Results and discussion of the global 
analysis for each reaction channel is presented in Sec. IV. 
Einally, conclusions are given in Sec. VH. 


II. REDUCED R-MATRIX THEORY 

A. The representations |54l . |55| for decay channel 
of 7n + l®0„ 

For the i^C(a, 7)1^0 reaction, the decay channel of 
7 „-|- 160 „ (n= 0 , I, 2 , 3 , 4 ) has the same total angular mo¬ 
mentum and parity J'"’ as the entrance channel i^C-|-a. 
Three angular momenta are involved at this stage for 
the state of compound nucleus, the spins of the target 
(or residual nucleus) and of the incident (or emitting) 
particle and the orbital angular momentum of the inci¬ 
dent (or emitting) particle. There are two conventional 
ways of combining these for the decay channel 7n-l-^60„. 
One is the channel spin scheme in which the vector sum 
of the spins of photon (1-^=1) and residual nucleus 
(Jf) is firstly formed giving the channel spin s 

I.y -)- Jf = S. ( 2 ) 

Channel spin s and the orbital angular momentum 1 of 
the photon are then combined to form the spin of the 
compound nucleus 

s + l = Ji. (3) 

Finally a state in the representation is labeled by the set 
of quantities {a (l 7 Jf ) s/JjM}. 

The alternative coupling scheme is called Devons and 
Goldfarb the ‘^-representation” [d^, which is used 
widely in the R-matrix theory mentioned above [s^, Is^ . 
Here the spin of photon is combined with the correspond¬ 
ing orbital angular momentum to form total angular mo¬ 
mentum for the photon states. It gives 

I 7 -t 1 = L. 


( 4 ) 
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Note:y^^,y^,,a = (I^=l, 
- 3' 


7040 


7030 


7o20 


7oio ■ 


Jy= 0)=0 


7o13 

7oi5 7oi2 

7o14 ■>'011 


e e" 

V 


7oi3 7oio 

7oi2 
I e e 

L 


£-representation 


channel spin scheme 


FIG. 1: The electric multiploes transition processes to 

ground state (0*^) described by channel spin scheme and 73- 
representation for system. 


The compound state’s spin is then given by 

L-t-Jf=Ji. (5) 


And a state of compound nucleus is labeled by the set of 
quantities {a (l7Jf ) LJf JjM} in the scheme. 

These two coupling schemes represent the same phys¬ 
ical situation in two different representations. In the 73- 
representation, parity conservation implies that 

TTi = ( 6 ) 

with L the multipolarity and the P mode (1 = electric, 0 
= magnetic) of the gamma ray, and tt^, tt/, n-y are the par¬ 
ity of the initial state Ji , final state Jf and the intrinsic 
parity of photon, respectively. Electric multiploes corre¬ 
spond to parity (-1)^ and I = Til, while magnetic mul¬ 
tiploes correspond to parity (-1)^+^ and I = L 
Based on the expression of Eq. we can deduce the 
parity conservation in the channel spin scheme 


7Ti — 7T j7T>y7Ti 


TrfTTy{—iy, I = L ±1 


(7) 


where tt; is parity of orbit angular momentum for the exit 
channels. 

For the comparison of two coupling schemes, Fig. [1] dis¬ 
plays the all possible transition processes to the ground- 
state (J'' = 0+) of system. The ground-state tran¬ 
sition of 73-representation, set of components has one 
value for each decay, which is less than the two com¬ 
ponents in the channel spin scheme of the Fig. [1] and 
the relevant photon reduced-width amplitude of this pro¬ 
cesses only exists one parameter in the collision ma¬ 
trix. Only when the channel capture is considered in 
73-representation, this parameter can be spli t into inter¬ 
nal and asymptotic channel contributions [^, . So the 

channel spin scheme provides more subsets to denote the 
ground-state transition and cascade transition, which is 
extremely beneficial to the interpretation of the observed 
experimental data more reliable. 


B. Wave function of compound nucleus 


For the ^^C(a,7)^®0 reaction, the transition of com¬ 
pound nucleus initial radioactive decay to ground 
state and four bound states are regarded as two body 
particle reaction channels, denoted by (n=0, 

1, 2, 3, 4), and the reduced masses of these channels 
are represented by relativistic energy. It is not neces¬ 
sary to consider how to decay to the ground state finally, 
so the problem of ‘particles are created or destroyed’ is 
avoided. The sum of integral cross sections of each re¬ 
action channel 7 n+^® 0 „ is equal to the cross section of 
production. 

Owing to the advantage of the channel spin scheme, 
all the channels of the system are represented as c = 
asl, where s is the channel spin, I is the relative angular 
momentum of the interacting particle of entrance or exit 
channels, and a identifies the interacting particle pair. 
The primary wave function 41 can be unfolded with dif¬ 
ferent exit channels ipc- Furthermore ipc can be expended 
with level wave functions X\^ which have different total 
angular momentum, parity and E\. Finally is ex¬ 
pressed with Eq. 2.6 (page 283) in Ref. [48l| : 


c 


E X\"fXc 
Ex-E 


n't 


( 8 ) 


So the total wave function for the initial state of can 
be expanded by the complete orthogonal set, the coeffi¬ 
cient of the expanded formula represent the probability of 
different reaction channel of all sorts of resonance energy 
state. Eq. ([8]) demonstrates that, if the primary gamma 
decay 7 n+^® 0 „ as the independent two body reaction 
channel, the set of level wave function Xx: of the theo¬ 
retical model contains all types of 7 - transition, whether 
direct decay to the ground state transition or cascade 
transition. Using the Eq. ([8]), one can obtain the funda¬ 
mental R-matrix relation, the collision matrix, the cross 
sections and so on. 

For the channel 7n+^®Ora, the electromagnetic inter¬ 
action is long range, therefore contributions to the colli¬ 
sion matrix for radiative capture reactions can come from 
large distances. Thus, in addition to the internal contri¬ 
bution to the collision matrix, there should also be chan¬ 
nel contributions [s^. However, for the chief ground state 
transition, owing to the large binding energy (7.16 MeV) 
of the with respect to the a-l-^^C threshold, its wave 
function decreases rapidly when the radius is larger than 
a certain value. So the internal contribution is strongly 
dominant, and the external part can be neglected. The 
results in Refs. [ 2 ^ prove to be reasonable and ef¬ 
fective for this approximation, that the external contri¬ 
bution accounts for lower than 3% at 0.3 MeV. For the 
cascade transitions, a parameter for the final state can be 
used to characterize the direct capture process of these 
transitions(see the parameters table). So a global fitting 
for whole system can be done using the standard 
R-matrix formulae of Ref. with a suitable channel 
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radius. 


C. Mathematical formalism of RAC code 

The practical formulae of our RAC code are introduced 
from the literatures [STl - f^ . On the R-matrix and the 
reaction cross sections, the codes are strictly compiled in 
accordance with the formulae of classic literatures (d^ . 
without any approximation. 

Explicitly, the i?-matrix, which represents all the in¬ 
ternal information concerning the structure of the com¬ 
pound system, is defined as 

N 

Xfi 


where Iasi reduced-width amplitude 

of entrance and exit channel, respectively. The matrix 
is defined by its inverse 

= ( 10 ) 

where Ex is the position of resonance level, Aa^ is the en¬ 
ergy shift, Ea^ is the total reduced channel width, which 
can represent the contribution of all un-considered chan¬ 
nels, such as the ^^C(a,a 2 )^^C in our fit. The additional 
quantity appearing in Eqs. (113 is 

N 

Aa// ^ ^ 'Ya's'I'l^ash (H) 

asl 

where Sx^^ is the shift factor calculated at the channel 
radius, and Bx/j, is the constant boundary parameter. 

The literature of of Lane and Thomas (Page 273) (d^ . 
gives the correction formula on the level width Tac and 
level shift Aac: 


Tac = 2Pc7L/^c, (12) 

P4R°,P,)-S°{1-R°X) 2 
=--7ac. (13) 

I 


where 


4 = (14) 

Here, A is the level of c reaction channel. Pc is the pen¬ 
etration factor and notation zero is the constant back¬ 
ground. These formulae are workable only based upon an 
approximation of single level. RAC is the multi-channel 
and multi-level i?-matrix formula without constant back¬ 
ground. When calculating the width and shift of some 
levels, the calculated values of R-Matrix with remaining 
levels are taken as the constant background of the level. 
The observed width can be related to the physical re¬ 
duced width amplitudes with a formula as follow: 

rf/ = rAc(^i + ^^7L^)^^- (15) 

The total width for a state A is then the sum 

= (16) 


With the relation between T-matrix and U-matrix of 
formula 


T;i's'l',asl = e^^--^5c.,s'l',asl - Ui,,,, ^^^sl. ( 17 ) 

for a reaction going through a —>■ a' the angle-integrated 
cross section is given as. 


(2J-f 1) 


U2 


sl' s'lj 


( 2 /i + 1 ) ( 2/2 + 1 ) 


'T'j 

^ a.'s'V ,asl 


(18) 


and Ii and I 2 are the projectile and target spins, respec¬ 
tively. It should be noted that the above equation does 
not hold for charged particle elastic scattering. 

For the corresponding differential cross section for¬ 
mula, a more rigorous calculation is involved in Ref. [dsj l 

= ( 2/1 + 1 )( 2/2 + 1 ) E 

(19) 

where Ac^'s'y',oisu is the amplitudes of the outgoing waves. 


Ack' s' y' ,0(sy{^Q^' ) — j [ Ca'{Oct'^Scn' s' y' jCusy 
Krv 


J MU' m' 


Several new quantities have been introduced in Eq. ((201) 
to define the angular dependence of the cross section. 
The term —Ca'{0a') represents the Coulomb amplitudes, 


while is the spheric harmonics function. 

The transverse character of electromagnetic wave re¬ 
quires that the projections of intrinsic spins of photon can 
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not be zero. For the transition with Jf = 0, when v' = 0, 
the 7 spin projection is zero. So when using the Eq. m 
to calculate the angular distribution of 7 decay, as long 
as ignoring the loop for z/' = 0 , the calculation will not 
include the contribution of 7 spin projection component 
zero. When this method is adopted, the angle-integrated 
cross section can be described effectively in our fit, but 
the corresponding differential cross section is not fitted 
accurately. So in the actual work, the longitudinal contri¬ 
bution of photon is employed to give a precise description 
for the all available data. 


D. Covariance statistic and error propagation law 

The uncertainty determination of the extrapolated S 
factor requires an error propagation of all relevant fit pa¬ 
rameters through the fit function taking into account the 
covariances. The theoretical formula about error propa¬ 
gation for our i?-matrix model fitting is as following: 


y - yo = D(P - Po), 

( 21 ) 

Dki = {dyk/dPi)o. 

( 22 ) 


Here y refers to vector of calculated values, D to sensi¬ 
tivity matrix, P to vector of i?-matrix parameters. Sub¬ 
script 0 means optimized original value, k and i stand for 
fitted data and i?-matrix parameter, respectively. The 
covariance matrix of parameter P is 

Vp = (D+V-^D)-^ (23) 

Here V refers to covariance matrix of the data to be 
fitted, and its inversion matrix can be expressed as fol¬ 
lowing: 


II 

1 

> 

/vr^ 

0 \ 


V 0 

1 -id 

> 


where Vi, ¥2 - ■■ Vk refer to the covariance matrixes of 
the sub-set data, which are independent with each other. 
The covariance matrix of calculated values is 

Vy = DVpD+. (25) 

The sensitivity matrix is quite useful in eliminating re¬ 
dundant fit parameters and in understanding which fit 
parameters are the most effective on the low-energy ex¬ 
trapolation of the S factor as discussed below. 

Formula adopted for optimizing with i?-matrix fitting 
is 

= (?7 — y)~'^V ~^(?7 — y) => minimum. (26) 

Here t] refers to the vector of experimental data, y refers 
to the vector of calculated values. Using covariance 


statistics and error propagation law, it enables us to get 
accurate expected values and standard deviation of S fac¬ 
tor. 

In addition the Peelle Pertinent Puzzle (PPP) was cor¬ 
rected by the method used in . RAC was used to pro¬ 
duce the accurate (error= 1 %) ®Li(n,a) and ^°B(n,a) 
cross sections, for International Evaluation of Neutron 
Cross section Standards [H [H]. And RAC was com¬ 
prehensively compared with R-Matrix code EDA and 
SAMMY of USA [H, [H^, the results were highly iden¬ 
tical when the same parameters are used. To verify the 
performance of the R-matrix code of a spectrum, 
we repeated the analysis of Ref. [ 2 ^ using their input 
data, and the same results were obtained. In a word, it 
is proved that the code RAC is reliable. 


III. EVALUATION AND FIT OF 
EXPERIMENTAL DATA 

A. The construction of reaction channel 

The observed position, width or life time of lev¬ 
els are displayed in Fig. [2] up to the 17.5 MeV energy 
region from Ref. iBflll . The cascade transition data in 
Ref. @ and Ref. [42j| reveal several 7 -ray cascade transi¬ 
tions from the 1“, 3“ and 2+ states at = 7.12, 11.49 
and 11.51 MeV, respectively, which are considered in our 
fit. Except the level = 2“ = 8.872 MeV), the 

other 31 levels contain a-|-^^C reaction channel. Based 
on the level scheme of the nucleus, the R-matrix 
analysis of the compound nucleus considers one par¬ 
ticle entrance channel ^^C-|-a, eight particle exit channels 
as shown in the Table UTl The i?-matrix calculations were 
performed for J'^ = O"*" (four real levels, one background 
level), J'^ = 1 “ (five real levels, one background level), 
= 2+ and 3“ (seven real levels, and one background 
level), = 4+ (four real levels, one background level), 
= 5“ (two levels, one background level) and J'^ = 


TABLE II: The reaction channels in our fit including Q- 
values, radii, the maximum of the orbital angular momen¬ 
tum, and the data, respectively, in the R-matrix calculation. 
For the capture reaction, Magnetic L-pole radiation is weaker 
than the corresponding electric L-pole radiation significantly, 
so the Magnetic transitions are not considered in our fit. 


Channel 

Q{MeV) 

R{fm) 

^max 

Data 

a + ^'^C 

0.000 

6.5 

6 


70 -b ^®Oo 

7.162 

6.5 

3 

AD, Sg.s. 

7i + ^®Oi 

1.113 

6.5 

1 

•S'e.os 

72 + 

1.032 

6.5 

1 

•S'e.is 

73 + 

0.245 

6.5 

1 

•S'6.92 

74 + 

0.045 

6.5 

1 

S 7.12 


-4.438 

6.5 

2 

AD, a 

p-bi®N 

-4.968 

6.5 

2 

AD 


“Here AD is the abbreviation of angular distribution. 
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FIG. 2: (Color online) Level scheme of the nucleus [6^. All states relevant for the analysis are indicated. 
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6 + (two real levels). The available data sets cover the 
energy from E^.m. = 0.9 MeV to E^.m. = 7.5 MeV in 
this fit, so the parameters above this region are fixed at 
values determined from the Ref. [b^. 

B. Evaluation of experimental data 

The basic principle of the data evaluation is that the 
database can reflect the information of nuclear structure 
and nuclear reaction accurately and objectively, no mat¬ 
ter which is to use the original data or the appropriate 
amendment. i?-matrix fitting requires the experimental 
data covering full energy region with complete energy 
points and continuous values, especially in the resonance 
peak area with the different types of data. Reliable ex¬ 
perimental data subset should satisfy the following re¬ 
quirements: In the resonance peak area, the sum of S 
factor in different reaction channel should be equal to the 
total S factor; The peak position of the different types 
of data should be consistent within the range of error; 
The principal value of different groups should be consis¬ 
tent within the range of uncertainty; The width data of 
resonance peaks are matched to the implied width in¬ 
formation of the other data; The integral value of the 
differential data should be equal to the corresponding in¬ 
tegral data; The integral data of different groups should 
span a broad energy range with a number of data points 
and have a good match with each other. 

According to the principle of maximum likelihood, a fit 
to a dataset with many types and large amount of points 
needs to meet the approximate statistical distribution, 
so the revisions of some dataset are reasonable. If one 
experimental point deviates from the expectations obvi¬ 
ously, such as the residual error larger than three times 
of uncertainty, the error of this point can be enlarged 
with the Letts’ criteria (Scr criteria); In the same type 
of data, if the difference of principal value is far greater 
than their uncertainties, the error of corresponding data 
should be amplified in the fitting; If the principal value in 
one group data deviates from the expected value wholly, 
the normalization to this dataset is needed in the fitting; 
If one high precision dataset is selected as the standard 
data in the evaluation, then some data with systematical 
deviation should be normalized to the standard data. 


C. Iterative fit 

The fits to the data of system are iterated, to 
achieve internal consistency. A file is a fixed record of 
the original data, which is to provide the original statis¬ 
tical error for the fit. Another file is a dynamic data file 
recording the evaluation process, which role is to provide 
the actually used data in fitting and is updated in the 
iterative process. In the file the original relative data 
values are replaced with the new normalized value, and 
the systematic error values are updated by the standard 



FIG. 3: The flow chart for the procedure of iterative fit. 

deviation (STD) of the new calculation. And the sta¬ 
tistical errors are renewed with the original one at the 
beginning, but some of them are corrected according to 
the Letts’ criteria. The ratio of the corresponding data 
in there two files is the new scaling factor or normaliza¬ 
tion coefficient. The scaling factor is adjustable in RAC, 
which is recorded in the parameter file together with the 
new R-matrix parameters. Fig.[3]shows the flow chart of 
R-matrix iterative fit procedure. 

With the continuity of iterative fit, the variation of 
scaling factor becomes smaller and smaller, and the prin¬ 
cipal values of relative experimental data are closer to 
their expectations. Similarly, the R-matrix parameters 
(RMP), fitted values and their standard deviations be¬ 
come more accurate. At last, all calculated values tend 
to very slight fluctuations, and the approaches the 
minimum. 


IV. RESULTS AND DISCUSSION 

The following subsections give details of the different 
reaction channels included in this analysis. Although 
they are described individually, the fits to the different 
reaction channel data sets have been performed simulta¬ 
neously. 


A. The reduced a-width amplitude for the bound 
states 

At energies of astrophysical interest, direct cross 
section measurement of capture reaction, such as 
^^B(p, 7 )^^C, is very difficult because of the Coulomb 
barrier, but it can be derived by the proton spectro- 














































scopic factor and asymptotic normalization coefficients 
(ANC) from the transfer reaction 

based on distorted wave Born approximation (DWBA) 
analysis The S factor of ^^C(a, 7)^®0 at astro- 

physical energies arises largely from the high-energy tails 
of subthreshold states 2^(Ex = 6.92 MeV) and IC^Ex = 
7.12 MeV) of but the properties of these states are 
only weakly constrained by cross-section measurements 
at higher energies. The cross section of transfer reac¬ 
tions d) and CLi, t) provides an alternative way for 
extracting the reduced a-widths for these states of 
In this fit, the of 1)” and 2^ bound states are fixed to 
the weighted average of two new measurements Bin, 
and the other subthreshold states, the ja of 0)'' and 3)" 
are adopted by the literature value of Ref. [3l| . While the 
j-y of the four states could vary within their uncertainties 
of literature . 


B. Total S factor 

The available total S factor of ^^C(a, 7)^®0 have been 
obtained in four independent experiments [rl-fT^. Fig. 0] 
illustrates the corresponding fitted values. In general, 
the fits are perfect where all the energy levels are ac¬ 
curate described. The measurement of Schiirmann et 
al. [71, Ig in inverse kinematics using the recoil mass sepa¬ 
rator ERNA allowed to collect data with high precision in 
a wide energy range, which would make a good restric¬ 
tion to the extrapolation of ground transition, cascade 
transitions and the total S factor. The data of Ref. [ 3 ] 
have not given definite numerical value of three narrow 
peaks 2 ^, 4+ and 0 ^, and the author’s personal com¬ 
munication considers that the relative numerical value is 
difficult to be determined, so the excitation energies and 
partial widths are fixed by including in the dataset of 
pseudo cross section points which were assigned by 50% 
errors around the resonance peaks. 

In the peak region of the I^ energy level {E^ = 9.58 
MeV), the evaluation of data is as follow. We can 
learn the Stot ~ 97.0 keV b from the measurement of 
Schiirmann et al. Q. The evaluated Sg.s ~ 76.0 is from 
six groups of Sg.s. Estimating the Se.os ~ 1-0 keV b, Sr a? 
« 7.0 keV b and S 7 ,i 2 ~ 20.0 keV b from Ref. [2^ and 
Ref. [l^ and making assumption for Sr. 13 ~ 1.0 keV b. 
In view of the above, the sum of the partial S factor is 105 
keV b , which is 1.08 times of the results Stot = 97 keV 
b from the measurement of Schiirmann et al. Q, that is 
to say, the experimental value of the total S-factor is sig¬ 
nificantly less than the sum of the experimental partial 
S factor. Relevant Ref. Q accounts for that the maxi¬ 
mal systemical error is 6.5%, so we choose 1.03 as the 
normalization coefficient of the data Q in the careful 
exploration, then the data become Stot = 100 keV b, 
which is lower than the sum of partial S factors. The 
systematical study shows that the Sr. 92 of Ref. [l^ has 
an increasing trend, and the S 7 ,i 2 of this paper has a 
decreasing trend. When taking the normalization coef- 



FIG. 4: (Color online) Results of the best R-matrix fit for 
the Stot data from Schurmann 2005: [ 3 ], Schiirmann 2011: [^, 
Flag 2012: [^, and Fujita 2013: [lol ]. together with the decom¬ 
position into different energy level contributions. For com¬ 
parison, the results of Kunz 2002: [^. Hammer 2005: [i^ . 
Schurmann 2012: [^, and Xu 2013: [6^ are shown in this 
figure, respectively. 


ficient of Sr .92 and 87,12 as 1.00 and 0.95, respectively, 
then the sum of partial S-factor is approximately equal 
to 100 keV b. Theresore we can get a satisfied dataset 
which has complete types and numerical self-consistency 
for the main resonance peak 1^. These constitute the 
skeleton of the whole database for the fits. 

Another skeleton of the dataset is the data on the peak 
region of 23 at Ec.m. = 4.358 MeV (See Fig.|4|). The data 
of Schurmann et al. Q is obtained by the adding of their 
components Sg.s, Sr.oSj Sr. 13 , Sr ,92 and S 7 . 12 , and it is 
consistent with the total S factor of Schurmann et al. [ 3 | 
very well. All kinds of the data of Ref. Q are used as the 
standard data, and the normalization coefficient is 1.03. 

Recently, the ^^C(a, 7)^®0 cross sections of Flag et 
al. 0 have been measured at four energy points, Ec.m. 
between 1.00 and 1.51 MeV, and the ElO and E20 com¬ 
ponents were derived with an accuracy comparable to the 
previous best data obtained with HPGe detectors. This 
data are first employed in the Stot ht, which have great 
influences on the S'tot(0.3 MeV). In Ref. [l^, total cross 
section measurements for Ec.m.= 2.4 and 1.5 MeV were 
performed at KUTL by using a tandem accelerator. And 
our fit results are relatively close to the principal values. 

In current research on S factor of Ref. [13, at higher 
energies, i.e., at Ec.m. > 2.8 MeV, resonance parameters 
taken from Ref. [60j were used in their R-matrix fit, in 
which the published data at hig h energy, such as a cap¬ 
ture measurements of Ref. were neglected. So the 
high-energy resonances from Ec.m. = 5 MeV to Ec.m. = 
6 MeV are overestimated apparently (please see the fit of 
Sg.s.). In addition, one should note that in the analysis 
of Ref. [13, |63] there is a clear disagreement at ener¬ 
gies around Ec.m. = 3 and 4 MeV, where the calculation 
underestimates total cross-section. The latest results of 
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FIG. 5: (Color online) Results of the best R-matrix fits 
(Black line) for the a spectra from Tang 2010: 
Azuma 1994: and Zhao 1993: [ 2 ^, together with the 

decomposition into p- (dashed line) and f- wave (dotted line) 
contributions. For comparison, the best fits of Tang 2010: [ 2 ^, 
Azuma 1994: [^. and Schurmann 2012: [iO] scaled by the 
corresponding coefficients are shown. 


Ref. [ 4 ^ are consistent with the available experimental 
data, but the high energy data are not analysed in a sim¬ 
ilar way of the R-matrix fit. In Ref. [^(NACREII) the 
total and partial S factors are analyzed with the poten¬ 
tial model, where the S factor at 2'^ {E^ = 11.52 MeV) 
is underestimated by the calculation. 


FIG. 6: (Color online) Results of the best R-matrix fits 
for the ground state transitions Sg.s. data from Brochard 
1972: f^, Ophel 1976: Kettner 1982: fp], Redder 

1987: 3, Ouellet 1996: Il3li . Kunz 2001: [l5|. Assungao 
2006: 3, Makii 2009: [18|| . Schurmann 2011: [^, Flag 
2012: [ 3 , together with the decomposition into different en¬ 
ergy level contributions. For comparison, the results of Brune 
1999: [p, Kunz 2002: [pi. H ammer 2005: Ip, Schurmann 
2012: 2^, Oulebsir 2012: (3i|. and Xu 2013: [6^ are shown. 


surement of Azmua et al [p most likely represents the 
currently closest approximation to the true a spectrum. 




C. a spectrum 


The shape of the low-energy part of the /3-delayed a 
spectrum of is very sensitive to the a -I- reduced 
width of the 1)" subthreshold state and 1^ state of 
in turn, which dominates the low-energy p-wave capture 
5'Eio(0.3MeV) of ^^C(q;,7o)^®Oo. In this energy region 
the reduced a widths are determined by the a spectra 
and the angular distributions of ^^C(q;, q:)^^C, which re¬ 
sults in the competition with each other in the fit. As 
shown in Ref. | 68 i | , there exists a limitation by the use of 
^^C(a,a)^^C data in Ref. for obtaining reliable val¬ 
ues of S_eio(0.3 MeV). So the a spectrum may help to 

give a better confirmation of the reduced a width ampli¬ 
tude of 1)" and I 2 ■ Included in this analysis are the three 


independent a spectra data of Refs. [25l - l28| . and the nor¬ 


malization for probability spectrum is used in the prac¬ 
tice to reduce the influence of systematical errors. Fig. [5] 
shows the fit to the normalized a spectrum together 
with the decomposition into p- and f-wave contributions, 
which suggests a significantly negative interference of the 
bound li state with the broad 1^ state that leads to a 
second peak at Ec.m. = 1.1 MeV and a minimum in the 
vicinity of 1.4 MeV. The dotted line denotes the con¬ 
tribution of 3“ state, which perfectly compensates this 
negative interference. The fit concluded that the inea- 


For the ground state transitions, the secondary data 
of Fio and F 20 multipoles were used in the previous R- 
matrix analysis independently . In general, these 

secondary data were obtained from the Legendre poly¬ 
nomials fit (page 510) to the angular distributions 
of ^^C(q;, 7 o)^®Oo measured at many discrete energies. 
Two methods of analysis (phase fixed or free) are often 
applied, however, the derived S factors S£;io and SB 20 
are significantly different for the same 7 angular distri¬ 
butions, see the Fig.12 of Ref. [l^, especially for the 
SB20- 

In our ht, the Sg,s.=S£;io+S £;20 is used for the ground 
state transition, which the proportion of S^io and SB 20 
are determined by the R-matrix ht to the relevant 7 angu¬ 
lar distributions. Fig. | 6 ] shows the ht to (a,7o) 
data of ten independent measurements and the calcula¬ 
tions of previous works. It is worth mentioning that af¬ 
ter the experiments by Dyer and Barnes , Kettner et 
al [ 2 ^ , and Redder et al [1^ , a weighted average value cr 
of 47 ± 3 nb at resonance peak of 1^ was used to derive a 
cross section at low energy, this data plays a vital role for 
the determination of Sg.s., and can be regarded as a crite¬ 
rion for normalizing the experimental data. Fven though 
the Sg,s. of Kettner et al [23| deviates systematically from 
the other data in this resonant region, it is dispensable 
since it is the only one that has the data points at above 
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FIG. 7: (Color online) The calculation of the Seio with the 
best R-matrix fit. For comparison, previous fits of Azuma 
1994: 1^, Ouellet 199 6: Jl^ . Brune 1999: J^, Gialanella 
2001: Kunz 2002: |6^ . Hammer 2005: Ej Schiirmann 

2012: [40| . Oulebsir 2012: |3l]]. and Xu 2013: |66| are shown in 
this figure. Data points shown are taken from Dyer 1974: [lll| . 
Redder 1987: [T^], Ouellet 1996: [T^, Roters 1999: Kunz 

2001: Gialanella 2001: (^. Assungao 2006 : [l6l| . Makii 

2009: (l8l| . Schiirmann 2011: and Flag 2012: 



FIG. 8: (Color online) The calculation of the Se20 with 
the best R-matrix fit. For comparison, previous fits of Ouel¬ 
let 1996: [T^, Brune 1999: [2^, Kunz 200 2: 1^ . Hammer 
2005: [3, Dufour 200 8: ^Il| . Oulebsir 2012: [3l|. Schiirmann 
2012: [40| . Sayre 2012: [4ll| and Xu 2013: are shown in this 

figure. Data points shown are taken from Redder 1987: [l3| . 
Ouellet 1996: [T^, Kunz 2001: [3, Assungao 2006: [I^, Makii 
2009: [l^ . Schiirmann 2011: [^T^nd Flag 2012: [^. 

3.0 MeV. It is noteworthy that if the normalization co¬ 
efficient is fixed by a factor of 0.87, the data of 1^ peak 
region is consistent with the other data, meanwhile the 
data above 3.0 MeV is well consistent with the Sg.s. of 
Schiirmann et al. 0. 

In higher e nergy re gion, there are five independent ex¬ 
periments pRl - l^ .l^ covering the 1;^ {E^ = 12.44 MeV) 
and IJ {Ex = 13.09 MeV). All the available measure- 



FIG. 9: Fits to the ^^C(a, 7o)^®Oo angular distributions of 
Ouellet 1996: [H at E^.m. = 1.362 (A), 1.370 (B), 1.382 (C), 
1.576 (D), 1.578 (E) and 1.580 (F) MeV. 



FIG. 10: Fits to the ^^C(q:, 7o)'^®Oo angular distributions of 
Ouellet 1996: O at E^.^. = 1.777 (A), 1.979 (B), 2.172 (C), 
2.383 (D), 2.390 (E) and 2.570 (F) MeV. 



FIG. 11: Fits to the ^^C(q:, 7o)'^®Oo angular distributions of 
Ouellet 1996: [3 at E^.^. = 2.590 (A), 2.835 (B), 2.982 (C) 
and 2.985 (D) MeV.] 
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FIG. 12: Fits to the ^^C(a, 7o)^^Oo angular distributions of 
Makii 2009: [3 at E^.m. = 1.225 (A), 1.470 (B), 1.467 (C), 
and 1.591 (D) MeV. 


FIG. 13: Fits to the ^^C(a, 7o)^®Oo angular distributions of 
Dyer 1974: [3 at E^.^. = 2.180 (A), 2.420 (B), 2.560 (C), 
and 2.831 (D) MeV. 


ments of the relative ratio of peak cross sections cr(l;^) 
/ct (14 ), are tabulated in table 3 of Ref. [31: in which the 
largest deviation of Ref. [i^l data from the other three 
data, lower about 20%, are evident. So in our fit, the 
data of Ref. [3 near IJ resonance are corrected to the 
data at the peak of 1;^ with a factor 0.81. Then the cross 
sections are found to be in a good agreement with the 
absolute data from Ref. [3 if normalization corrections 
(maximum of ±20%) are applied. The remaining data 
from Ref. [3> HH13 show good agreement in the shape 
of the excitation curves, and the normalization factors 
are given in Table. IIIII 

The corresponding R-matrix calculation of angular dis¬ 
tributions for the reaction ^^C(a, 7 o)^®Oo are illustrated 
in Fig. |9]- Fig. [22] at representative a energies from E^.m. 
= 1.002 to Ec.m. = 6.075 MeV. The angular distribution 
measured at the 2^ (E^ = 9.84 MeV), has the familiar 
E20 pattern, is symmetric with respect to 90°, while the 
distributions obtained at other energies are asymmetric 
about 90°, clearly indicating the presence of both ElO 
and E20 amplitudes in the capture mechanism. With 
the much improved 7 ray angular distributions in our R- 
matrix calculation, it can now be possible to derive more 
accurate values for the cross sections of the Eio and E 20 
transitions to the ground state of Fig. [Tjand Fig. |5| 
show the calculations of the S eio and S ^20 together with 
all the available experimental data. Although the data 
of S^jio and 8520 are not used in the fits, these data lie 
in two sides of our calculation uniformly, which in turn 
illustrates the rationality and self-consistency in our R- 
matrix fit of the angular distributions and Sg.s, ground 
state transitions. 

In Refs. [3 [3 tfie subthreshold state and the 
I 2 resonance may interfere destructively and result in a 
significantly lower Sbiq. However, the constructive solu¬ 
tion is strongly favored and the destructive interference 
pattern has been eliminated in our calculation of angu¬ 
lar distribution, resulting in a value of Seio(0.3 MeV) = 



angle 9 (degree) 

FIG. 14: Fits to the 7o)^®Oo angular distributions of 

Redder 1987: [3 at E^.^. = 1.710 (A), 2.360 (B), 2.684 (G), 
and 2.830 (D) MeV. 


98.0 ± 7.0 keVb. 


The cross section around the Ec.m. = 2.5-3.0 MeV 
region is a rapidly changing function of energy, which 
strongly depends upon the interference scheme between 
the resonance 2^ and other E 20 amplitudes. But the rel¬ 
ative E 20 -E 20 interference sign is not well determined by 
the integral capture data, i.e., the best result of Ref. [3 
in an interference pattern determined by the high-energy 
data of Ref. Q is different from most previous analyses 
Refs, [illnl- The interference scheme has been com- 
mendably constrained by the angular distributions calcu¬ 
lation in our R-matrix fitting near this resonance, which 
are in accordance with the new measurement result of 


Ref. [3: and the extrapolation values of S£; 2 o( 0.3 MeV) 
= 56.0 ± 4.1 keVb in our calculation. 
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FIG. 15: Fits to the ^^C(a, 7o)^^Oo angular distributions of 
Assungao 2006: 0 at Ec.m. = 1-310 (A), 1.340 (B), 2.268 
(C), 2.660 (D), 2.677 (E) and 2.684 (F) MeV. 


FIG. 18: Fits to the 7o)^®Oo angular distributions 

of Fey 2004: [l^ at Ec.m. = 2.343 (A), 2.455 (B), 2.578 (C), 
2.607 (D), 2.652 (E) and 2.682 (F) MeV. 
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FIG. 16: Fits to the ^^G(a, 7 o)^^Oo angular distributions of 
Flag 2012: [| at Ec.m. = 1-002 (A), 1.308 (B), 1.416 (C) and 
1.510 (D) MeV. 
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FIG. 17: Fits to the 7 o)^®Oo angular distributions 

of Fey 2004: [13 at E^.m. = 1.666 (A), 1.965 (B), 2.040 (C), 
2.116 (D), 2.192 (E) and 2.230 (F) MeV. 


FIG. 19: Fits to the ^^C(a, 7 o)'^®Oo angular distributions of 
Fey 2004: [13 at E^.m. = 2.684 (A), 2.687 (B), 2.757 (G) and 
2.780 (D) MeV. 
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FIG. 20: Fits to the 7 o)'^®Oo angular distributions of 

Kunz 1997: [3 at Ec.m. = 2.250 (A), 2.400 (B), 2.480 (C) 
and 2.684 (D) MeV. 
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FIG. 21: Fits to the ^^C(a, 7o)^^Oo angular distributions of 
Larson 1964: [l^ and Kernel 1971: [2l| at Ec.m. = 5.295 (A), 
5.565 (B), 5.910 (C), 6.000 (D), 5.775 (E) and 5.910 (F) MeV. 
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FIG. 22: (Color online) Fits to the ^^C(q, 7o)^®Oo differen¬ 
tial cross-section data at 45° (A), 61° (B), 90° (C), and 135° 
(D) from Larson 1964: [T^, Kernel 1971: [^ . Ophel 1976: 
and Mitchell 1964: [^. 


E. 


FIG. 23: (Color online) Results of best R-matrix fit for 

the cascade transitions Se.os from Schiinnann 2011: @] and 
Matei 2006: [IJ , together with the decomposition into differ¬ 
ent energy level contributions. For comparison, the results of 
Schiirmann 2011: [^, Matei 2006: [Ij, and Xu 2013: [6^ are 
shown in this figure. 


cited state (07 6.05 MeV). But Ref. Q concludes that 
the Se.os is negligible in the energy region less than 3.3 
MeV, so it only gives the experimental data above the 
energy. 

In our fit, the data of Ref. Q is regarded as standard 
data and the normalization coefficient is 1.03. The energy 
regions of the data in Ref. [2^ and Ref. Q have overlap 
around 3.5 MeV. The data of Ref. can be normal¬ 
ized by that of Ref. Q, and the normalization factor is 
0 . 88 . This forms a dataset of Se.os which covers full en¬ 
ergy region with complete energy points and continuous 
values. This transition can therefore be estimated to be 
Se.os (O.SMeV) = 4.9 ±1.2 keV b, where S^;! is the most 
important contribution in Se.os (0.3 MeV). The Se.os ob¬ 
tained by this work is from the systematic analysis of the 
whole system. Hence compared with previous anal¬ 
ysis, our result is much firmly based on the experiments 
and is reliable. 


Radiative a-particle capture into the first excited, = 
0+ state at 6.05 MeV excitation energy have been inves¬ 
tigated recently in two independent experiments sm, 
however, there exists big differences of Se.os (0.3 MeV) 
each other. So it is necessary to give a detailed re¬ 
search and discussion. In the work of C. Matei et.al. (^ . 
86 . 05 ( 0.3 MeV) = 25ti6 b is obtained by fitting the 
experimental results therein. It mainly comes from S^ji 
and partially from Se 2 and is with large error. In contrast 
to the analysis of Ref. [ 2 ^ , the extrapolation of Ref. @ 
suggests a negligible contribution from this amplitude, 
85 . 05 ( 0.3 MeV) < 1 keV b by analyzing their data, which 
is mainly contributed by Se 2 while little by S^i. Ref. 
and Ref. Q use the same experimental method, and both 
their original data show the 7 contribution of the first ex- 


F. ^2)1602 

Very little data exists about the transition into the 
= 6.13 MeV state (J’^ = 3“) except for the 2;^ resonance 
at Ex = 11.60 MeV and 3^ resonance at Ex = 11.52 
MeV of Ref. Q . The parameters of these resonance can 
be sufficient to describe this data and the fit result in 
86 . 13 ( 0.3 MeV) = 0.2 ± 0.1 keVb. 

G. i2C(a,73)i®03 

Four cascade data of 86. 92 [E El, El, Ell cover a range 
from Ec.m. = 1-4 to 5.5 MeV, and the cross section 
at astrophysical energy is largely governed by the di- 
































14 



FIG. 24: (Color online) Results of the best R-matrix fit 
for the cascade transitions Se.is from Schiirmann 2011: @] 
together with the decomposition into different energy level 
contributions. For comparison, the result of Xu 2013: is 

shown in this figure. 



FIG. 25: (Color online) Results of the best R-matrix fit for 
the cascade transitions S ^.^2 from Redder 1987: Kettner 

1982: [^, Kunz 2001: [l^, and Schiirmann 2011: [g, together 
with the decomposition into different energy level contribu¬ 
tions. For comparison, the results of Schiirmann 2012: 
and Xu 2013: are shown in this figure. 


reel capture process, from s-, d- and g-wave captures in 
Refs. la. fl^ 1^. With the reasonable normalization 
of Refs. a good fit to these experimental data are 

achieved from resonance parameters and a direct capture 
parameter for =2^ (see parameter table), resulting in 
86 . 92 ( 0.3 MeV) = 3.0 ± 0.4 keV b, which is consistent 
with the result of Ref. [i^. Fig. [25] shows the results of 
85.92 as well as its decomposition into different level con¬ 
tributions. Also, the normalization factors are given in 
Table HH 



FIG. 26: (Color online) Results of the best R-matrix fit for 
the cascade transitions S'r.iz from Redder 1987: [T^], Kunz 
2001: [ 1 ^, and Schiirmann 2011: [§], together with the decom¬ 
position into different energy level contributions. For compar¬ 
ison, the results of Schiirmann 2012: and Xu 2013: 

are shown in this figure. 

H. i2C(a,74)i®04 

The capture of 87,12 would be expected to proceed 
mainly via p- and f-wave direct process and resonance 
transition at low energies [1, [T^. With resonance pa¬ 
rameters and a direct capture parameter for J’^=l)”, a 
good fit to the experimental data [H, [T^ is obtained 
(see parameter table). And the extrapolated S factor for 
this transition is also small, 87 , 12 ( 0.3 MeV) = 0.6 ± 0.2 
keVb. The normalization factors of these applied data 
are given in Table HIT] Fig. [55] shows the results of 87,12 
as well as its decomposition into different level contribu¬ 
tions. 


I. ^^C(a,ay^C 

The most widely used elastic scattering data of a- 
particles on contain the a particle information for 
all relevant states in which can be obtained with 
rather high accuracy. Previous elastic-scattering data 
have been used to determine the scattering phase shifts 
for individual angular momenta Refs. [H, and so on. 
Such a procedure is necessary in cases when the analysis 
of ^^C(a, 7)^®0 is restricted to only one particular an¬ 
gular momentum, but the interference structures in all 
data, associated with all the resonance states, have been 
neglected |4^ . Taking all angular momenta into ac¬ 
count simultaneously, R-matrix fits of four group angu¬ 
lar distributions and the associated data (32 - 1^ are pre¬ 
sented from Ec.m. =1.1 to Ec.m. = 5.85 MeV in Fig. 1571 
to Fig. [38] with the strictly theoretical formulae of Ea.ITQl 
In order to reduce the space of the paper, the hgures from 
Ec.m. = 5.85 MeV to Ec.m. = 7.5 MeV of Refs. [H [s^ 
are not shown in this paper. 
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FIG. 27: Fits to the ^^C(a,a)^^C angular distributions of 
Plaga 1987: [H at E^.m. = 1.100 (A), 1.480 (B), 1.550 (C), 
1.555 (D), 1.704 (E), and 1.870 (F) MeV. 


FIG. 30: Fits to the angular distributions of 

Plaga 1987: [H at E^,m. = 2.829 (A), 2.888 (B), 3.038 (G), 
3.131 (D), 3.196 (E), and 3.284 (F) MeV. 
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FIG. 28: Fits to the ^^G(a,f:v)^^G angular distributions of 
Plaga 1987: at E^.m. = 2.003 (A), 2.153 (B), 2.228 (G), 

2.302 (D), 2.303 (E), and 2.378 (F) MeV. 


FIG. 31: Fits to the ^^C(o',f:v)^^G angular distributions of 
Plaga 1987: [H at E^.m. = 3.338 (A), 3.488 (B), 3.638 (G), 
3.788 (D), 3.944 (E), and 4.088 (F) MeV. 
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FIG. 29: Fits to the ^^C(a,a)^^C angular distributions of 
Plaga 1987: at E^.m. = 2.452 (A), 2.457 (B), 2.528 (C), 

2.588 (D), 2.678 (E), and 2.738 (F) MeV. 


FIG. 32: Fits to the angular distributions of 

Plaga 1987: [H at E^.^. = 4.125 (A), 4.163 (B), 4.200 (G), 
4.238 (D), 4.275 (E), and 4.290 (F) MeV. 
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FIG. 33: Fits to the ^^C(q,q:)^^C angular distributions of 
Plaga 1987; [H at E^.m. = 4.352 (A), 4.364 (B), 4.387 (C), 
4.424 (D), 4.462 (E), and 4.499 (F) MeV. 



angle 0 (degree) 

FIG. 35: Fits to the ^^0(0;, a')^^G angular distributions of 
Plaga 1987: [H at E^.^. = 4.769 (A), 4.844 (B) and 4.919 
(G) MeV. 




FIG. 34: Fits to the ^^G(a,a)^^G angular distributions of 
Plaga 1987; [H at E^.m. = 4.507 (A), 4.542 (B), 4.582 (G), 
4.619 (D), 4.657 (E), and 4.694 (F) MeV. 


FIG. 36: Fits to the q:)^^G angular distributions of 

Tischhauser 2009: [H at E^.m. = 2.291 (A), 3.192 (B), 3.913 
(G) and 4.902 (D) MeV. 


The scattering data by Plaga et al. [s^ were obtained 
in the considerably better energy range in comparing 
with the other studies. Differential cross-section data for 
all 35 angles in the range 9iab = 22° to 163° and for 51 
energies from E^.m. = 1.0 to 4.9 MeV are included in the 
fit. Data points in the vicinity of narrow resonances are 
also contained for this analysis. Level parameters from 
the fit are in an excellent agreement with those reported 
in Ref. . The interference structures in the data, asso¬ 
ciated with resonance states in the energy range covered 
by this data, are well reproduced by the R-matrix fits. 
All the results are shown in Figs 113 to Figsl35] 

Recently the angular distributions of ^^C(a,a)^^C in 
the a-energy range of 2.6-8.2 MeV, at angles from 24° 
to 166° have been measured at the University of Notre 
Dame using an array of 32 silicon detectors [H, 1^ . 
The relative differential cross-section excitation curves 
for eight selected detector angles and the four angular 
distributions for energies near the E^.m. = 2.291 (1^), 
3.192 (4j), 3.913 (2^), and 4.902 (Oj) MeV resonances 


are available. To reduce the amount of computations, 
only the four angular distributions are employed in this 
fit, which the angular distributions are found to be in 
a good agreement with those data. Fits are shown in 

FiglMl 

The best quality a-scattering cross-section data above 
proton separation energies are shown in FigITfland Figl551 
of Ref. |35l| , which have good coherence to experimental 
data. 


J. i=C(q,qi) 12C and i=C(a,p)^®N 

Good fits of ^^C(a, q:i)^^C and are help¬ 

ful to reduce the uncertainty produced by the distant lev¬ 
els, and then to improve the fit precision of ^^C(a, 7)^®0 
S factor subsequently. Previously, the only available an¬ 
gular distribution data of the ^^C(q:, ai)^^C reaction has 
been obtained at incident energies Ec.m. = 5.963, 6.015 
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FIG. 37: Fits to the ^^C(q,q:)^^C angular distributions of 
Morris 1968: [H at E^.m. = 4.950 (A), 5.163 (B), 5.224 (C), 
5.243 (D), 5.298 (E) and 5.320 (F) MeV. 



FIG. 38: Fits to the ^^G(a, a')^^G angular distributions of 
Morris 1968 [H] at E^.m. = 5.375 (A), 5.550 (B), 5.625 (G), 
5.775 (D), 5.825 (E) and 5.850 (F) MeV. 



0 45 90 135 180 

angle 0 (degree) 


FIG. 40: Fits to the ^^G(a,ai)^^G angular distributions of 
Mitchell 1965: [13 at E^.m. = 5.963 (A), 6.015 (B) and 6.105 
(G) MeV. 



E^___(MeV) 



FIG. 39: (Golor online) Fits to the ^^G(q,qi}^^G cross sec¬ 
tion of Mitchell 1965: js^ and deBoer 2012: [s^ . 


FIG. 41: Fits to the ^^0(0:, ai)^^G differential cross section 
data of Mitchell 1965: [13. 

and 6.105 MeV in Ref. [13. And the excitation curves of 
^^C(a,ai)^^C and ^^C(q;,p)^^N for four selected detec¬ 
tor angles have been measured in the same energy range 
in the paper. But an absolute scaling was not reported 
in these measurements. In Ref. (ssl. [^. new yield-ratio 
data for the reactions ^^C(a,ai)^^C and ^^C(q;,p)^®N 
were performed at the University of Notre Dame in or¬ 
der to provide additional data for a comprehensive R- 
matrix analysis of compound-nucleus reactions populat¬ 
ing The data are in the form of yield ratios where the 
^^C(a,ao)^^C yields measured at an angle 61iab = 58.9° 
are used as the reference data. And the transformational 
angular distribution data and cross sections are obtained 
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by a private communication with Dr. R. J. deBoer. 

The R-matrix fits of ^^C(a, ai)^^C angle-integrated 
cross section are shown by the solid lines in Figl39l il¬ 
lustrating the consistency level of the simultaneous fit to 
the data of Ref. together with the relative data in 
Ref. [S^- Fits of the available angular distribution data 
are shown in Fig 001 to Fig 031 in which the significant 
contributions originate from states at Ex = 12.95 (2^), 
13.13 (Sg ), and 13.27 ( 3lD MeV. Fits for the angular dis¬ 
tribution data of Ref. for the reaction ^^C(q;,p)^®N 
are shown in Fig ]44l and the normalization factors are 
given in Table Hill 


FIG. 42: Fits to the ^^C(a,ai)^^C angular distributions of 
deBoer 2012: [H at E^.m. = 5.964 (A), 5.984 (B), 6.001 (C), 
6.020 (D), 6.039 (E) and 6.059 (F) MeV. 



angle 0 (degree) 

FIG. 43: Fits to the ^^G(a,ai)^^G angular distributions of 
deBoer 2012: HI at E^.m. = 6.096 MeV (A), 6.113 (B), 6.133 
(C), 6.153 (D) and 6.169 (E) MeV. 


V. SUMMARY 


This study presents a new R-matrix theory for the 
^^C(a,7)^®0 S factor at helium burning temperature, 
and a number of applications to demonstrate the appli¬ 
cability and versatility of this theory. The final result 
of S(0.3 MeV) = 162.7 ± 7.3 represents the most pre¬ 
cise extrapolation of the ^^C(a,7)^®0 S factor at helium 
burning temperature based on a set of complementary 
data including all available information of system so 
far. This is to our knowledge the first published analysis 
meeting the precision requirements on ^^C(a, 7)^®0. The 
whole S factor from 0.3 MeV to 10 MeV provides astro- 
physical reaction rate of ^^C(a, 7)^®0 with a sound basis 
for researches of nucleosynthesis and evolution of stars. 



5.8 6.0 6.2 5.8 6.0 6.2 

Ec.m. (MeV) 


FIG. 44: Fits to the ^^C(q:,p)^®N differential cross section 
data of Mitchell 1965: [37| . 
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TABLE III: Scaling factors for datasets which have no reported absolute scale, and the normalization coefficient for the absolute 
data. 
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TABLE IV: Parameters obtained from the R-matrix fits, particle partial widths for the levels, and the radiation widths for the 
7 -ray transitions considered in this work compared with those from literatures. 


Res. no. 

E^MeV) 

TTT 

Ex{MeV) 

a(7i + I 2 ) 

(s,0 

7c..i(MeV)i/" 

r(keV) 

Ref. [60] 

01 

6.049 

ot 

-0.1113 

a+^'^C 

(0,0) 

-1.37(18) X 10“" 

Fixed 

19.7(55) 





7o+^®Oo 

(1,1) 

2.14(41) X 10”® 

6.88(50) X 10-^2 

6.86(50) X 10-^2 

02 

12.049 

ot 

4.865 
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4.99(151) X 10“® 

1.02(25) X 10° 
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Fixed 
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5.91 X 10“^ 

Fixed 

Ttot = 166(30) 






(1,1) 

-2.85 X 10“^ 
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05 

BG°- 

ot 
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a-h^^C 

(0,0) 

2.23(9) X 10° 







7o-l-^®Oo 
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-3.01(20) X 10“'‘ 
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7.119 

ir 

-0.4515 

a+i^C 

(0,1) 

-8.08 X 10"^ 

Fixed 

62.0(170) 





7o-l-^®Oo 

(1,0) 

2.20(107) X lO""* 

5.7(3) X 10"® 

5.5(3) X 10"® 





7o+^®Oo 

(1,2) 

1.61(46) X 10”^ 







7i+^®Oi 

(1,0) 

-2.06(898) X 10“° 

3.0(1) X 10"“ 

< 3.0 X 10"®° 





72 + ^®02 

(2,1) 

-7.32(83) X lO-'^ 

4.5(1) X 10"® 

4.6(10) X 10"® 





73 + ^®03 

(1,0) 

5.54(1300) X 10"° 

5.5(2) X 10"“ 

< 1.0 X 10"° 





74 + ^®04'’ 

(0,1) 

-1.00(67) X 10-2 



07 
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^2 

2.295 

a-h^^C 
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3.37(3) X 10"^ 

322(5) 
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-7.37(25) X 10"® 

7.25(23) X 10"° 

1.56(12) X 10"® 





7o-|-^®Oo 

(1,2) 

-2.99(9) X 10"° 







7i+'®Oi 

(1,0) 

-1.39(15) X IQ-'^ 







73 + ^®03 

(1,0) 

-6.46^5) X 10"® 

5.8(8) X 10"'^ 

1.4(14) X 10"° 





74 + ^®04 

(0,1) 

1.36(68) X 10"® 

1.5(15) X 10"° 

7.8^6) X 10"° 

08 

12.442 

la" 

5.289 

a+i^C 

(0,1) 

1.13(2) X 10"^ 

162(7) 

102(4) 





7o+^®Oo 

(1,0) 

-8.15(137) X 10-'‘ 

2.3(3) X 10"® 

1.2(2) X 10-2 





7o+^^Oo 

(1,2) 

-2.73(13) X 10-2 







7i+^®Oi 

(1,0) 

-5.25(59) X 10-'‘ 

1.2(3) X 10-'‘ 

1.2(6) X 10"® 





72-|-^®02 

(2,1) 

1.39(51) X 10"® 

3.4(25) X 10"® 

7.0(30) X 10"®" 





74-|-^®04 

(0,1) 

2.89(55) X 10"® 

9.0(34) X 10"® 

1.3(5) X 10"®" 





ai+i^C 

(2,1) 

1.12(4) X 10"^ 

2.9(2) X 10-2 

2.5 X 10-2 





p-b^®N 

(1,0) 

1.16(10) X 10"^ 

1.8(3) X 10"^ 

9.0(1) X 10"® 

09 

13.088 

ir 

5.859 

a+i^C 

(0,1) 

-4.71(49) X 10-2 

27(4) 

45(18) 





7o+^®Oo 

(1,0) 

-3.56(52) X 10"® 

1.20(4) X 10-2 

3.2(5) X 10-2 





7o+^®Oo 

(1,2) 

-2.42^4) X 10-2 







7i+^®Oi 

(1,0) 

-7.77(118) X IQ-'^ 

2.7(8) X 10-“^ 

2.4(5) X 10"® 





72 + ^®02 

(2,1) 

-3.54(86) X 10"® 

2.8(13) X 10-'‘ 

4.0(20) X 10"®" 





74 + ^®04 

(0,1) 

-9.53(136) X 10"® 

1.3(4) X 10"® 

1.35(40) X 10"® 





ai+i^C 

(2,1) 

8.30(52) X 10-2 

8.5(9) X 10"^ 

1.0 X 10° 





p-b^®N 

(1,0) 

2.35^3) X 10"^ 

3.4(4) X 10^ 

1.1(2) X 102 

10 

17.510 

15 

9.710 

a+i^C 

(0,1) 

-8.99(387) X 10-2 

Fixed 

29(9) 

11 

BG 

le" 

11.434 

a-b^^C 

(0,1) 

1.23(1) X 10-2 







7o+^®Oo 

(1,0) 

-5.75(65) X 10"^ 







7o+^®Oo 

(1,2) 

7.40(237) X 10"® 







7i+^®Oi 

(1,0) 

-9.71(133) X 10"^ 







p-b^®N 

(1,0) 

-8.95(290) X 10-2 



12 

6.917 


-0.2449 

a+i^C 

(0,2) 

1.68 X 10"® 

Fixed 

26.7(103) 





7o-|-^®Oo 

(1,1) 

1.99(8) X 10"® 

9.9(3) X 10"® 

9.7(3) X 10"® 





7o+^®Oo 

(1,3) 

2.81(37) X 10"^ 







7i+^®Oi 

(1,1) 

-6.75(5) X 10-* 

2.7(3) X 10"® 

2.7(3) X 10"® 





72 + ^®02 

(2,0) 

1.30(29) X 10"® 

9.0(30) X 10"° 

9.0 X 10"° 





73 + ^®03‘' 

(1,1) 

1.96(23) X 10-2 




“BG is the abbreviation of background level. 
^Direct capture for the S7.12. 

^Ref. [si . 

^Direct capture for the 56.92- 







21 


TABLE V: Parameters obtained from the R-matrix fits, particle partial widths for the levels, and the radiation widths for the 
7-ray transitions considered in this work compared with those from literatures. 


Res. no. 

El^MeV) 

TTT 

E;A(MeV) 

a(7i + I 2 ) 

(s,Z) 

7c.si(MeV)i/" 

r(keV) 

Ref. [60] 

13 

9.844 

2+ 

2.684 


(0,2) 

1.34(2) X 10"^ 

0.71(19) 

0.625(100) 





7o+^®Oo 

(1,1) 

-1.62(3) X 10"^ 

2.0(1) X 10-° 

5.7(6) X 10-° 





7o+^®Oo 

(1,3) 

-1.03(3) X 10"= 







7i+'®Oi 

(1,1) 

-3.27(43) X lO-'^ 

4.4(12) X 10-'^ 

1.9(4) X 10-° 





73 + ^®03 

(1,1) 

-1.03(9) X 10"® 

2.1(4) X 10-° 

2.2(4) X 10-° 

14 

11.520 

2+ 

4.314 

a+^^C 

(0,2) 

6.86(3) X 10"^ 

74(1) 

71(5) 





7o+^®Oo 

(1,1) 

1.77(17) X 10“® 

6.8(2) X 10-° 

6.1(2) X 10-° 





7o+^®Oo 

(1,3) 

-2.02(15) X 10“^ 







7i+^®Oi 

(1,1) 

1.64(13) X 10“® 

3.1(5) X 10-° 

3.0(5) X 10-° 





72 + ^®02 

(2,0) 

2.46(16) X 10“^ 

2.2(3) X 10-® 

2.0 X 10-®“ 





73 + ^®03 

(1,1) 

-1.66(14) X 10"^ 

2.0(3) X 10-° 

2.9(7) X 10-° 





74 + ^®04 

(2,0) 

-1.90(32) X 10"'‘ 

1.1(4) X 10-® 

< 5.0 X 10-° 

15 

13.020 

2+ 

5.833 


(0,2) 

8.15(20) X 10“^ 

112(5) 

150(10) 





7o+^®Oo 

(1,1) 

-8.51(66) X lO"'^ 

2.2(2) X 10-° 

7.0(20) X 10-° 





7o+^®Oo 

(1,3) 

-7.54(74) X 10“^ 







ai+^^C 

(2,0) 

-4.04(37) X 10"^ 

4.1(8) X 10-° 

5.0(20) X 10-°“ 






(1,1) 

-8.23(45) X 10“^ 

1.6(2) X 10° 

1.5(2) X 10°“ 

16 

15.90 

2+ 

8.300 


(0,2) 

2.22(9) X 10”^ 

Fixed 

Ttot = 600 





7o+^®Oo 

(1,1) 

-1.34(8) X 10"® 

Fixed 

Pcr^/rTOT = o.4ey 

17 

16.443 

2^ 

9.281 


(0,2) 

-3.37(26) X 10"^ 

Fixed 

Pa/PTOT = 0.28 





7o+^®Oo 

(1,1) 

1.38(7) X 10"® 

Fixed 

PcP-y/ProT = 0.45ey 





ai+^^C 

(2,0) 

-9.50 X 10"^ 

Fixed 

Ptot = 22(3) 






(1,1) 

8.75 X 10“® 

Fixed 


18 

17.129 

2| 

9.967 

a+i=C 

(0,2) 

-7.68 X 10-2 

Fixed 

P^/Ptot = 0.37 





ai+^"C 

(2,0) 

6.62 X 10“® 

Fixed 

Tror = 107(14) 






(1,1) 

1.88 X 10-2 

Fixed 


19 

BG 

2+ 

22.618 

Q-b^^C 

(0,2) 

2.22(2) X 10° 







7o+^®Oo 

(1,1) 

1.15(14) X 10-° 







7o+^®Oo 

(1,3) 

1.67(48) X 10-° 







7i+^®Oi 

(1,1) 

1.01(12) X 10-2 







73 + ^®03 

(1,1) 

3.18(21) X 10-° 







ai+^'C 

(2,0) 

-4.90(42) X 10° 








(1,1) 

-1.35(58) X 10° 



20 

6.1299 

3r 

-0.1032 

Q-b^^C 

(0,3) 

-7.39(17) X 10-2 

5.44(25) 

2.35(80) 





7o+^®Oo 

(1,2) 

5.83(31) X 10-° 

2.47(7) X 10-® 

2.60(13) X 10-® 

21 

11.600 

^2 

3.969 

Q+l^C 

(0,3) 

3.14(2) X 10-° 

718(10) 

800(100) 





72 + ^®02 

(2,1) 

5.38(3930) X 10-° 

2.4(346) X 10-® 

1.0 X 10-®“ 





73 + ^®03 

(3,0) 

-1.05(59) X 10-° 

2.8(31) X 10-° 

1.0 X 10-°“ 





74 + ^®04 

(2,1) 

6.13(540) X 10-° 

1.6(29) X 10-° 

2.0 X 10-®“ 

22 

13.129 

33“ 

6.082 


(0,3) 

1.02(3) X 10-° 

69(5) 

90(14) 





7o+^®Oo 

(1,2) 

-1.63(79) X 10-° 

9.1(30) X 10-° 

1.0 X 10-° 





72 + ^®02 

(2,1) 

-1.83(30) X 10-2 

7.5(25) X 10-° 

8.0 X 10°“ 





ai+^'C 

(2,1) 

-2.88(30) X 10-° 

4.26(32) X 10° 

2.09(6) X 10° 






(1,2) 

-1.12(119) X 10-° 

9.8(6) X 10-° 

1.0 X 10° 

23 

13.259 

34“ 

6.063 

Q-b^"C 

(0,3) 

-3.77(21) X 10-2 

13(4) 

9(4) 





72 + ^®02 

(2,1) 

-1.79(15) X 10-2 

8.0(13) X 10-° 

9.2(15) X 10-°“ 





ai+^'C 

(2,1) 

3.44(15) X 10-2 

2.96(36) X 102 

8.2(11) X 10° 






(1,2) 

4.43(12) X 10-° 

1.6(1) X 10° 

4.1 X 10° 

24 

14.100 

35 

7.149 

Q-b^^C 

(0,3) 

-5.154 X 10-° 

Fixed 

150(75) 





ai+^^C 

(2,1) 

-1.03(63) X 10-2 

Fixed 

Ptot = 750(200) 

25 

15.408 

36 

8.246 

a-b^^C 

(0,3) 

-2.57 X 10-’^ 

Fixed 

T^/Ttot = 0.58 





ai-b^'C 

(2,1) 

2.24 X 10-° 

Fixed 

Ptot = 133(7) 






(1,2) 

9.83 X 10-° 




“Ref [H]. 
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TABLE VI: Parameters obtained from the R-matrix fits, particle partial widths for the levels, and the radiation widths for the 
7 -ray transitions considered in this work compared with those from the literature. 


Res. no. 

E^^MeV) 

j: 

£;A(MeV) 

a(7i + I 2 ) 

{s,l) 

7,.i(MeV)i/2 

r(keV) 

Ref. 

26 

15.828 

3^ 

8.266 

a+^'^C 

(0,3) 

-2.24 X 10“^ 

Fixed 

Fa/FTOT = 0.21 





ai+^"C 

( 2 , 1 ) 

1.03 X 10“^ 

Fixed 

Ftot = 703(113) 

27 

BG 

38 

22.618 

a-b^^G 

(0,3) 

-1.64(1) X 10° 







74+^®04 

( 2 , 1 ) 

-1.58(63) X 10"^ 



28 

10.361 

4t 

3.196 

Q-b^"C 

(0,4) 

2.17(3) X 10"^ 

28.4(6) 

26(3) 





7o+^®Oo 

(1,3) 

1.15(2650) X 10"'‘ 

5.6(20) X 10-“ 

5.6(20) X 10-“ 





72 + ^®02 

(4,0) 

4.27(47) X 10“® 

4.7(10) X 10-'^ 

< 1.0 X 10-° 





73 + ^®03 

(3,1) 

-4.25(33) X 10”° 

5.2(8) X 10-° 

6.2((6) X 10-° 

29 

11.094 

4 

3.934 

Q+l^C 

(0,4) 

9.78(108) X 10-° 

0.296(32) 

0.28(5) 





72 + ^®02 

(4,0) 

1.05(19) X 10“^ 

3.7(1.3) X 10-° 

3.1(13) X 10-° 





73 + ^®03 

(3,1) 

7.04(82) X lO-'^ 

2.7(6) X 10-° 

2.5(6) X 10-° 

30 

13.879 


6.835 

Q-b^"C 

(0,4) 

-5.00(56) X 10-2 

50.3(49) 

Fq/Ptot = 0.65(5) 





ai+^"c 

( 2 , 1 ) 

-2.01(14) X 10-^ 

Ftot = 74(7) 

Ftot = 77(7) 

31 

16.844 

‘it 

9.682 

Q-b^"C 

(0,4) 

1.06(21) X 10-^ 

Fixed 

Ftot = 567(60), Fc/Ptot = 0.28 

32 

BG 

it 

30.230 


(0,4) 

2.10(3) X 10° 







73 + ^®03 

(3,1) 

-3.44(30) X 10-2 







ai+^"C 

( 2 , 2 ) 

-8.50(118) X 10-^ 



33 

14.660 

sr 

7.498 

a+^^C 

(0,5) 

3.47 X 10-^ 

Fixed 

Ftot = 672(11), Fc/Ftot = 0.94 

34 

16.900 

^2 

9.748 

a+^'^C 
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2.9 X 10-^ 

Fixed 

= 700 
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Ftot = 70(8), Tc/Ttot = 0.28 
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9.133 
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( 0 , 6 ) 

2.83 X 10-^ 
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Ttot = 422(14) 
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